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1.0  INTRODUCTION 


Historical  records  have  been  kept  for  World  Meteorological  Organization 
stations  worldwide  for  a  number  of  years.  Because  of  their  volume,  these 
records  have  been  summarized  In  various  forms  to  be  readily  available  to  a 
user.  One  method  of  summarization,  or  compaction.  Is  by  means  of  an 
empirical  cumulative  distribution  function  (cdf ) .  The  empirical  cumulative 
distribution  function  is  simply  the  tabulated  cumulative  relative  frequencies, 
or  probabilities  that  a  given  climatological  variable  will  fall  below 
specified  values.  Somerville  and  Bean  (1979)  have  shown  that  a  number  of 
climatological  variables  may  be  modeled  with  closed  form  distribution  functions. 
For  a  given  location  and  time  (e.g.,  month  and  hour)  the  historical  observations 
can  be  used  to  obtain  estimates  of  the  model  parameters. 

In  the  present  study,  we  have  attempted  to  model  visibility  for  locations 
in  Germany  for  which  there  are  no  historical  records.  This  has  been  done  by 
first  developing  models  for  stations  with  existing  records,  and  then  attempting 
to  express  the  parameters  of  the  models  as  a  function  of  such  measured 
variables  as  elevation,  relative  elevation,  proximity  to  water,  etc.  using 
regression  techniques. 

2.0  THE  USE  OF  THE  WEIBULL  DISTRIBUTION  FOR  VISIBILITY  PROBABILITIES 
The  Weibull  distribution  has  been  found  to  be  useful  in  modeling 
visibility  by  Somerville,  Bean  and  Falls  (1979).  The  probability  that 
visibility  is  less  than  x  miles,  using  the  Weibull  distribution  is  given  by 
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0 

F(x)  -  1  -  e-ax 

where  a  and  0  are  parameters.  The  values  for  the  parameters  a  and  0  vary 
depending  on  location,  time  of  day  and  month.  For  West  Germany,  the 
parameters  were  estimated  by  Somerville  and  Bean  (1981)  for  eight  three- 
hour  periods  during  the  day  and  for  each  month  for  30  locations.  The 
object  was  to  use  the  estimates  of  a  and  0  to  determine  an  empirical 
relationship  of  the  parameters  with  other  measurable  variables  using 
regression  techniques.  The  stations  used  in  the  study  are  listed  in  the 
Appendix. 


3.0  METHODS  FOR  OBTAINING  VISIBILITY  MODELS  FOR  LOCATIONS  WITHOUT 

HISTORICAL  RECORDS 

There  are  many  variables  which  may  influence  visibility.  Variables 
investigated  in  this  study  include  the  following. 

1.  Elevation. 

2.  Relative  elevation.  This  is  the  elevation  divided  by  the  "average 
elevation"  of  the  surrounding  area.  The  "average  elevation"  is  calculated 
by  forming  a  circle  of  radius  20  kilometers  about  the  location  in  question 
and  obtaining  the  elevations  for  a  grid  of  values  on  the  circle. 

3.  East-west  elevation  difference.  This  is  the  difference  between  the 
elevation  20  kilometers  east  and  20  kilometers  west  of  a  station. 

4.  North-south  elevation  difference.  This  is  the  difference  between 
the  elevation  20  kilometers  north  and  20  kilometers  south  of  a  station. 


3 


5.  Proximity  to  major  bodies  of  water.  This  is  an  indicator  variable 
which  has  the  value  1  if  the  location  is  near  (within  30  kilometers)  to  an 
ocean,  large  river  or  lake,  and  has  the  value  0  otherwise. 

6.  Population  density  (individuals  per  100  square  kilometers). 

7.  Relative  humidity  (for  a  given  month  and  hour  period). 

8.  Mean  wind  speed  (for  a  given  month  and  hour  period). 

9.  Mean  precipitation  (by  month). 

10.  Latitude. 

11.  Longitude. 

12.  Functions  of,  and  interactions  between  all  of  the  above. 

For  a  given  month  and  hour  period,  the  values  for  a  and  8  given  in 
Somerville  and  Bean  (1981)  for  30  stations  in  Germany  were  regressed  on  the 
values  for  the  above  variables.  That  is,  each  of  a  and  B  were  modeled  as  a 
function  of  variables  thought  to  be  candidates  for  estimating  visibility 
probabilities.  Using  these  models,  values  for  a  and  6  can  be  estimated 
and  probabilities  for  visibilities  obtained  provided  the  variables  in  the 
regression  equations  are  available. 

A  stepwise  regression  program  was  used  to  ascertain  which  of  the  above 
variables  could  be  used  as  predictors  for  a  specified  month  and  hour  period. 
Stepwise  regression  is  a  model  building  procedure  which  sequentially  includes 
"statistically  significant"  variables  in  a  model,  entering  the  most  important 
variables  first,  and  removing  previously  included  variables  which  appear  to 
be  redundant  given  new  variables  have  been  Included,  in  an  iterative  fashion. 
For  a  more  complete  description,  the  reader  is  referred  to  Draper  and  Smith 
(1981). 
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For  each  month,  and  for  each  of  the  eight  3-hour  periods  expressions 
for  each  of  a  and  8  were  obtained  which  were  functions  of  a  selected  (by 
the  stepwise  regression  program)  set  of  the  above  variables.  The  regression 
equations  were  then  re-determined,  using  the  selected  variables  in  a 
specially  designed  non-linear  regression  program,  which,  using  the  data  from 
all  30  stations  simultaneously,  selected  the  regression  coefficients  in  the 
expressions  for  a  and  8  which  minimized  sum  of  squares  of  the  differences 
between  the  observed  and  predicted  probabilities  over  all  distances  and  all 
stations.  In  non-linear  regression,  it  is  necessary  to  have  reasonable 
"starting  values"  for  the  coefficients  to  be  determined,  which  are  then 
refined  by  iteration.  The  "starting  values"  were  the  coefficients  in  the 
original  regression  equations  for  a  and  8  • 

The  expression  which  was  minimized  was 


30  14  r  -1 

Z  Z  E,(x.)  -  F.(x, ;  a, ,8.) 
j=l  i=l  L  3  1  3  1  3  2 


(3.1) 


r  th 


Here  x^  is  the  i  visibility  distance  included  in  the  RUSSWO's  (the  14 
distances,  in  miles,  were  1/4,  5/16,  1/2,  5/8,  3/4,  1,  5/4,  3/2,  2,  5/2,  3, 
4,  5  and  6).  The  index  j  corresponds  to  the  station.  (x^)  is  the 
empirical  probability  that  the  visibility  is  less  than  x^  miles  for  the 
j**1  station  as  obtained  from  the  RUSSWO. 


F^xjcyBj)  =  1  -  e~ajX  j 

is  the  model  cumulative  distribution  function  (Weibull)  for  the  j 


(3.2) 

th  station 


the  eight  3-hour  periods  are  0000-0200,  0300-0500,  etc.,  and  are  referred 
to  as  3-hour  periods  1,  2,  3,  ...  .  All  times  are  LST. 
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where 


6.  =  d,  +  6.2  +  ...  +  6  z  . 

j  1  1  ij  m  mj 


(3.4) 


The  y  ,  ...»  y^  are  the  values  at  the  j1"*1  station  of  the  k  variables 

selected  by  the  stepwise  regression  program  to  obtain  a,  and  z^,  z  ,  ...  z^ 
are  the  values  at  the  station  of  the  m  variables  selected  by  the  stepwise 

regression  program  to  obtain  B.  The  values  y^,  y^ . yfc  and  d^,  5  ,  ...,  d^ 

are  the  "improved"  regression  coefficients  for  the  regressions  of  a  and  3  on 
the  selected  variables. 

Exhibit  3.1  illustrates  the  steps  in  the  model  building  procedure  for 


a  given  month  and  three-hour  time  period  for  the  "Variables  Model." 


2. 


Using  the  Weibull  model,  where 

8 

F^xjOj.gj)  «  1  -  e~ajX  3 
choose  constants  Yq.  y^,  ...,  y^  and  6q,  6^,  ..., 

i  2 


l  l  |E  (xi 

i  i  L 


)  -  »jUli.1.8j) 


' 


6  so  as  to  minimize 
m 


where  the  indices  i  and  j  refer  to  the  i1"*1  distance  and  the  jttl  station, 


and  Ej (x)  is  the  empirical  cumulative  distributive  function  for  j 


th 


station. 

In  the  iterative  non-linear  regression,  use  the  variables  selected  in 
Step  2.  Initial  values  for  the  y's  and  6's  are  those  obtained  in  Step  2. 


Exhibit  3.1 

MODEL  BUILDING  STEPS  FOR  A  SPECIFIED  MONTH  AND  HOUR  PERIOD 


FOR  "VARIABLES  MODEL" 
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4.0  RESULTS 

For  the  30  German  stations  used  in  the  study,  the  values  of  a  and  8, 
for  each  of  the  96  combinations  of  month  and  hour  period,  are  given  in 
Somerville  and  Bean  (1981). 

For  stations  for  which  no  historic  records  are  available,  but  for 
which  values  or  estimates  may  be  used  for  such  predictor  variables  as 
average  relative  humidity,  average  windspeed,  elevation,  etc.,  Exhibit  4.3 
lists  the  formulas  for  obtaining  values  of  a  and  6  to  be  used  in  the  models. 

The  variables  used  in  the  formulas  are  functions  of  those  given  in 
Exhibit  4.1 


ABBREVIATION 
MP 
EL 
AE 

NE 
SE 
EE 
WE 
RH 

WA 

WN 

Exhibit  4.1 


ORIGINAL  VARIABLE  NAME 

Mean  precipitation  in  inches  for  a  given  month 
Location  elevation  in  feet 

Average  elevation  in  feet  (for  locations  20  kilometers 
distant 

Elevation  20  kilometers  north  of  the  location 

Elevation  20  kilometers  south  of  the  location 

Elevation  20  kilometers  east  of  the  location 

Elevation  20  kilometers  west  of  the  location 

Mean  relative  humidity  for  a  given  month  and  hour  period 
(values  between  0  and  100) 

Indicator  variable  for  proximity  to  water 
(1  if  major  water  body,  0  otherwise 

Mean  windspeed  for  a  given  month  in  knots 


ABBREVIATIONS  FOR  ORIGINAL  VARIABLE  NAME 
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SYMBOL  USED  IN  FORMULAS 

TRANSFORMATION  ON  ORIGINAL  VARIABLE 

RH 

RH3/105 

EL 

el3/io9 

MP 

mp3/iooo 

WN 

wn3/io 

WA 

WA 

DN 

(NE-SE)2/105 

DE 

(WE-EE)2/105 

ER 

(EL/AE)2 

Exhibit  4.2 

CODED  VARIABLE  USED  IN  FORMULAS 


Originally,  the  stepwise  regression  program  included  the  variables  in 
Exhibit  4.1  plus  their  squares,  cubes,  reciprocals,  reciprocals  of  their  squares 
and  two  and  three  way  interactions  of  those  terms.  Experimentation  indicated 
that  in  most  cases  if  some  power  of  a  variable  appeared  in  a  prediction 
equation,  the  other  powers  contributed  little.  Also,  in  many  cases,  one  power 
appeared  about  as  useful  as  another  for  prediction  purposes  (e.g.,  square  vs. 
cube) .  Thus  on  the  basis  of  many  runs  a  suitable  power  was  chosen — or  equiva¬ 
lently  a  variable  transformation  was  made  so  that  the  resulting  prediction 
equations  would  appear  more  elementary.  For  example,  use  of  the  cube  of  the  value 
for  relative  humidity,  and  the  square  of  the  ratio  of  elevation  to  average 
elevation  were  apparent  simplifications.  In  addition,  to  be  able  to  limit  the 
range  of  the  values  for  the  coefficients,  and  to  have  some  commonality  of  ranges, 
the  variables  were  coded. 

As  an  example,  RH  in  all  of  the  formulas  in  Exhibit  4.3  actually  refers 
to  10  multiplied  by  the  cube  of  mean  relative  humidity  corresponding  to  a 
specific  month  and  hour  period.  Exhibit  4.1  gives  the  original  variable  names 
and  their  abbreviations.  Exhibit  4.2  gives  the  meaning  of  the  symbols  used  in 


the  formulas  in  Exhibit  4.3. 
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M 

0  H 
N  0 


T 

H 

U 

R 

DEP 

VAR 

EQUATION 

1 

1 

ALPHA 

0.957<SD-01  +  0 . 21340+01 

* 

EL*  HP 

+ 

0,54600-03 

* 

WN*WA 

BETA 

0 . 12640+01  +  -0 . 82290-03 

* 

UN*DN 

RMS 

0.6271019E-01 

1 

2 

ALPHA 

-0.43380+00  +  0 i 84560-01 

* 

RH 

+ 

0.20590-03 

* 

WNtON 

BETA 

0  , 11690+01  +-0. 25650-01 

* 

DN*ER 

RMS 

0 , 814438PE-01 

1 

3 

ALPHA 

-0.40220+00  +  0.84410-01 

* 

RH 

+ 

0,50730-03 

* 

F.L*WN 

BETA 
+  -0 . 

0,11060+01  +-0. 27920-01 
37030-03  *  WNtON 

* 

ELfER 

+ 

0,81150-03 

* 

WN*UA 

RMS 

0 , 7493435E-01 

1 

4 

ALPHA 

0.17530+00  +  0.50470-03 

* 

FL*MN 

+ 

0. 12380+00 

* 

EL*WA 

BETA 

+-0. 

0,11760+01  +-0, 25790-02 
62450-02  *  UN*WA 

* 

RH*DN 

+ 

-0,12970-02 

* 

El.*WN 

RMS 

0 . 7366557E-01 

1 

5 

ALPHA 
+  0. 

-0.99820-01  +  0 , 4000D-01 
41300-07  *  MN*UA 

* 

RH 

+ 

0.14470-03 

* 

El.tWN 

BETA 
+  -0. 

0.12210+01  +-0. 41780-01 
53110-03  *  WN*DN 

* 

RH*WA 

+ 

-0.46300-02 

* 

FI.*OE 

RMS 

0 , 5516620E-01 

1 

6 

ALPHA 

-0.15410+00  +  0.48950-01 

* 

RH 

+ 

0,35640-02 

* 

UN*UA 

BETA 

+  -0  < 

0.12020+01  +-0. 36020-01 
31250-03  *  UN*DN 

* 

RH*WA 

+ 

-0.48990-02 

* 

F.l.*OE 

RMS 

0 , 591071 IE-01 

1 

7 

ALPHA 

0.13800+00 

BETA 
+  -0. 

0.11550+01  +  0,73290-01 
26890-01  *  DN+ER 

* 

EL  *FR 

+ 

0,31130-02 

* 

WN*MA 

RMS 

0 . 9058049E-01 

1 

8 

ALPHA 

0,87840-01  +  0,15740-01 

* 

EL 

BETA 

0.14430+01  +-0. 18670+00 

* 

ER 

+ 

-0.54680-03 

* 

WN*ON 

RMS  0  ,  6419713E-01 


Exhibit  4.3 

EQUATIONS  FOR  a  AND  B  (VARIABLES  MODEL) 


M 

0  H 

w  n 

t  II 

H  ft 


MAR  RQIJA1  ION 

ALPHA  0 « 74950-01  +  0 
+  0 » 37880-02  *  WA*DE 


.  13420-01  *  FL*ER  +-0. 23100+00  *  MPtOE 


BETA  0.12630+01  +-0, 84910-03  *  UNtDN 

RMS  0  .5034752E-01 

2  ALPHA  0.96730-01  +  0,45170-03  *  EL*WN 

BETA  0,12350+01  +-0. 30580-01  *  ON*ER 

RMS  0.8766603F-01 

3  ALPHA  0. 11830+00  +  0,27790-03  *  RH*UN 

BETA  0 ♦ 1 0470+01  +-0, 46930-03  *  WNCDN 

RMS  0 , 4498680E-0 1 


ALPHA  0,12140+00  +  0.24370-03  #  RH*WN  +  0.34640-01  *  WA*FR 
BETA  0 » t023D+01 


RMS  0.5699237F-01 

5  ALPHA  0.51550-01  +  0, 
+  0 , 73190-04  *  UN*ON 

BETA  0.13510+01  +-0, 
+-0,75460-03  *  WN*DN 

RMS  0 , 3508572E-01 

A  ALPHA  0.40050-01  +  0, 
+  0,41740-04  *  WN40N 

BETA  0 . 13AA0+01  +-0. 
+-0,44080-03  *  WN*DN 

RMS  0.372599AF-01 

7  ALPHA  -0.83080-01  +  0. 

+  0.10500-03  #  WN*OM 

BETA  0.14600+01  +-0. 
+  0,58110-02  *  RH*ER 

RMS  0. 4414B44E-01 

8  ALPHA  0.30660-01  +  0, 

+  0.48520-02  *  WA*DE 


.13800-03  *  RHtUN  +  0. 
+  0,13530-03  *  WN*OE 


13400-01  *  RH*WA 


63490-01  *  RH*WA  +-0. 11050-01  *  FL *BE 


.14320-03  *  RHtUN  +  0. 
+  0.11180-03  *  WN*DE 


11020-01  *  RH*WA 


.62320-01  *  RH*WA  +-0. 88890-0?  *  FL.&OE 


.26360-01  *  RH  +  0. 
+  0.50580-01  *  WA*ER 

.61190-01  *  RHtWA  +-0 . 
+-0, 48720-03  *  FLtUN 


29760-04  *  FL*MN 


64540-02  *  RHtON 


.26320-02  *  RKtEL  +  0.93760-04  *  RHtUN 


BETA  0.13670+01  +  0.38190+00  *  FL*MP  +-0. 76720-03  *  WK*ON 
RMS  0 , 5248782F-01 


Exhibit  4.3  (Continued) 


t  0  oep 

H  R  VAR  EQUATION 

3  1  ALPHA  0,14470-01  +  0.139211  +  03  *  F_L*MP 

BETA  0,14400+01  +  0,13800  +  0.3  *  MPfcWA 

RMS  0.4449725E-01 

3  2  ALPHA  0,47490-01  +  0. 59040+00  *  RH*HP 

BETA  0,13920  +  01  +-0. 35420-01  *  FI.*FR 

RMS  0 , 4278142E-01 

3  3  ALPHA  0.92140-01  +  0.42530-01  *  E!.  *WA  +  0.35040+01  *  MPtER 

BETA  0,11910+01  +  0,33580-02  *  UN  +-0. 12430-02  *  EL*UN 

+-0 . 1325D+00  *  EL*WA  +-0. 48450-03  *  WN#ON 

RMS  0 , 3995951E-01 

3  4  ALPHA  0.72910-01  +  0.27080-01  *  Fl.tWA  +  0.84200-03  »  F.l.*OE 

BETA  0,13710+01  +  0,54470-02  *  RH*WA  +-0. 14500-02  *  EL*UN 

+-0. 35000-03  *  UH+ON  +  0.10310-02  *  WHtER 

RMS  0 , 3992903E-01 

3  5  ALPHA  0.20900-01  +  0.78700-03  *  RH*F.I.  +  0.59790-02  *  RH*WA 

BETA  0,14340101  +-0. 20880-01  *  RH*WA  +-0.t3340-02  *  FI.KUN 

RMS  0 . 3074222E-01 

3  4  ALPHA  0.10990-01  +  0.13840-02  *  RH*Fl.  +  0. 31840-03  *  WN*WA 

BETA  0,17580+01  +-0,45720-01  *  El.*WA  +-0. 77820-01  *  EL*EF: 

RMS  0 . 2724449E-01 

3  7  ALPHA  0.15050-01  +  0. 29280+00  *  El.*KP  +  0.30010-01  *  El  *WA 

BETA  0 . 1714D+01  +  0,83480-01  *  RH*ER  +-0. 40040-02  *  F.LKUN 

+-0. 11410+03  *  MPtWA  +-0. 23470-01  *  DN*FR  +  0.33130-01  *  OF*F.R 

RMS  0 , 4043101E-01 

3  8  ALPHA  0.14710-01  +-0.1447D+01  *  HP  +  0.10730+01  *  EI.CMP 

BETA  0,194/0+01  +-0,20420-02  *  F.LHWN  +-0. 55430+00  *  EL*UA 

+-0, 45910+01  *  MP*DN  +-0, 52080-01  *  WA*OE  +  0.87010+00  *  WA*ER 

RMS  0 , 337394tE-01 


s 

I 


12 


M 

0  H 

N  0 

T  U  OEP 

H  R  VAR  EQUATION 

4  1  ALPHA  0.56600-02  +  0.76670-03  *  RHTEL 

BETA  0.21390+01  +  0,38790-02  *  DE  +-0. 75320-02  *  ELTUN 

RMS  0.3333507E-01 

4  2  ALPHA  -0 . 50350-01  +  0,15730-01  *  RH 

BETA  0,14960+01  +-0,22970-02  *  EI.TWN  +  0,72030-01  *  HATER 

RMS  0 .58953 18E-01 

4  3  ALPHA  -0 • 3018D-01  +  0.17420-01  #  RH 

BETA  0,17910+01  f-0.8281.D-0t  *  RH  +-0,73890-03  »  EI.TWN 
+-0,48640-03  *  WNTDN 

RMS  0 . 4748061E-01 

4  4  ALPHA  0.22530-01  +  0.13500-01  *  EL*WA  +  0.46350-03  *  EL  TOE 

BETA  0,13460+01  +  0,10440+00  *  RH  +-0. 18160-02  *  ELTUN 

+-0. 43680-03  *  WNTDN 

RMS  0, 304941 IE-01 

4  5  ALPHA  0.60490-02  +  0,49480-03  *  RH*EL  +  0,52240-02  *  El.TUA 

BETA  0,19360+01  +-0. 37730-01  *  El.TER 

RMS  0 . 2361841E-01 

4  6  ALPHA  0.17730-02  +  0,38230-04  *  RHTWN  +-0. 29070-03  *  RHTPN 

+  0,17370-02  *  El.TUA  +  0,19260  +  00  *  NPTER  +  0.55290-03  *  DNTER 

BETA  0.20030  +  01  +-0, 62730-01  *  ELTER 

RMS  0 , 1719512E-01 

4  7  ALPHA  0.16040-02  +  0.14680-01  »  MPTWN  +  0,98040-03  T  PNTER 

BETA  0,21870+01  +  0,13050-01  *  RHTWA  +-0.22H0+00  T  MPTUN 

+-0,38490-02  *  UNTDN 

RMS  0. 2852185E-01 

4  8  ALPHA  0,51500-02  +  0.39160-04  T  ELTON 

BETA  0,21650+01  f-0, 11290-02  *  WNTDN  +-0. 12890-02  »  WNTOE 

RMS  0 . 3007792E-01 


Exhibit  4.3  (Continued) 
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M 

0  H 

N  0 

T  U  OEP 

H  R  VAR  EQUATION 

5  t  ALPHA  0 . 74100-02 

BETA  0 . 1 8970401 

RMS  0.34t3559E-01 

5  2  ALPHA  0.42420-01 

BETA  0 » 13240401  4  0,21370-01  *  RHtWA 

RMS  0 .  4455042E-01 

5  3  ALPHA  -0  •  15550-01  4  0 , 1  0450-01  *  RH 

t  BETA  0,20080401  4-0,1.0321)400  *  RH  4-0.70110-03  *  ELTUN 

4-0  «  54890-03  *  WNTON 

RMS  0 . 4337838E-01 

5  4  ALPHA  0 . 45240-02  4  0.44940-04  *  EL*UN  4  0.11980-01  *  EL*WA 

BETA  0  >  2302D40 1  4-0.42170-02  *  EL*WN  4-0.34140400  *  EL*WA 

4-0,44070-03  *  UNION 

RMS  0.2932565E-01 

5  3  ALPHA  0 . 2198D-0?  4  0.31370-05  *  El.TWN 

BETA  0,21380401  4  0,14250-02  *  El. TUN  4-0.47450-03  *  UNTON 

RMS  0 . 1477727E-01 

5  4  ALPHA  0.87550-03  4  0.18950-04  *  El.tDN  4  0.21430-02  *  MPtUN 

BETA  0.23180401  4-0.19340400  *  MP*WN 

RMS  0. 1045017F-01 

5  7  ALPHA  0.44920-03  4-0.27950-05  *  RH*DN  4  0.44270-04  t  ELtDE 

4  0.29580-02  *  MPKHN  4  0.4912D-05  *  WN*DN  4-0.10550-03  *  WA*DN 

BETA  0 . 2467D401  4  0.15180-01  *  EI.*WA  4-0.27420400  t  MPtUN 
4-0,14990-02  *  WNTDN 

RMS  0 . 1484205E-01 

5  8  ALPHA  0,34430-02 

BETA  0,22040401  4-0.25140-02  *  WNKOE 

RMS  0.2483483E-01 


Exhibit  4.3  (Continued) 


14 


u 

OFP 

R 

OAR 

EQUATION 

1 

ALPHA 

0 . AA09D-02 

+  0.87870-04 

* 

ER 

BETA 

0.20090+01 

+  0.44190+00 

* 

KP*DF 

RMS 

0 , 3904  4A0F- 

•01 

2 

ALPHA 

0.50630-01 

BETA 

0.12770+01 

+  -0. 38370-03 

* 

FL  *WN 

RMS 

0 .4901203E- 

-01 

3 

ALPHA 

0.39680-01 

+  0,20980-01 

* 

F.L*WA 

BETA 

0,98480+00 

+  0.82330-01 

* 

RH 

RMS 

0 » 4973479F- 

01 

4 

ALPHA 

0.23540-02 

+-0,33270-03 

* 

UN 

+  0,25780+00  *  MP*WA  +  0.77A9D-04  * 

FETA  0.20850+01  +-0,43590-01  *  RH*ER 

RMS  0.2040592F-01 

ALPHA  0 . 91 240-03  +-0,10550-03  *  UN 
+  0 » 45200-01  *  NP*WA 


♦-0. 42230-03  *  WN*DN 


+  0 . 1 9A2D-03  *  RH*WN 
DN*ER 


+  0,81770-04  t  RH*UN 


BETA  0.23810+01  +-0,63580-01  t  RHtF.R 


RMS  0.12504B3E-01 

A  ALPHA  -0.11480-02  +  0,14240-0?  *  RH  +  -0 . 56970-04  *  ON 

+  0 » 33540-05  *  WN*0N  +  0.32850-04  *  WA*DN 

BETA  0 .23120+01  +  0,48200-02  *  RHtER 

RMS  0 » 12A5939F-01 

< 

7  ALPHA  0 » 12AAP-02  +-0.8088D-04  *  ON  +  0,61700-04  *  EL*WN 

+  0.30750-05  *  WN*DN 

BETA  0.22320+01  +-0. 41570-02  *  EI.*WN 

RMS  0 . 2325A18F-01 

8  ALPHA  0.21810-02  ♦  0.57080-04  *  RH*F.L  +  0.47240-04  *  EL*WN 

BETA  0.23250+01  +-0. 51850-02  *  EL*UN 

RMS  0.2444602E-01 


Exhibit  4.3  (Continued) 
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M 

0  H 
N  0 

T  U  OEP 

H  R  VAR  EQUATION 


7  1  ALPHA  0,48940-02 

BETA  0 , 1 9740+01 

RMS  0.4442354E-01 


7 

2 

ALPHA 

-0.45320-01  +  0 i 3  2880-01 

* 

RH 

BETA 

0,13870+01 

RMS 

0.7033592E-01 

7 

3 

ALPHA 

-0,17410-01  +  0 . 99460-02 

* 

RH 

+  0.14870-01 

t 

EL.tUA 

BETA 

0 , 16540+01  +-0 «  36440-0 1 

* 

RH 

RMS 

0.4588298E-01 

7 

4 

ALPHA 

0,88320-03  +  0.10410-03 

* 

RHBUN 

BETA 

0. 25360+01  +-0 <  79800-02 

* 

RH*WN 

RMS 

0.2934230E-01 

7 

5 

ALPHA 

0.98430-03  +-0. 57080-04 

* 

UN 

+  0. 53760-04 

* 

RH*UN 

BETA 

0  » 23570+01  +  0,28340+00 

* 

WA 

+-0. 77850-01 

* 

RHBER 

RMS 

0 , 1 29081 3E-01 

7 

6 

ALPHA 

+-0,: 

0 , 1 6630-02  +-0 i 62030-03  * 
28290-03  *  El.BER  +  0,14160-< 

El.  BMP  +  0.14140-04 
>2  *  MPBWN 

* 

EL  +  UN 

BETA 

0.19340+01  +  0,25040-02 

* 

EL*WN 

RMS 

0 , 8789920E-02 

7 

7 

ALPHA 

0.11090-02  +-0, 34200-03 

* 

WH 

+  0.18940-03 

« 

RH*WN 

BETA 

0 , 19840  101  +-0, 11020-02 

* 

ELBWN 

RMS 

0. 1794515E-01 

7 

a 

ALPHA 

0 , 1 7240-02 

BETA  0.25260+01 
RMS  0 « 3249431F.-01 


16 


M 

0  H 
N  0 


T 

U 

DEP 

H 

R 

VAR 

EQUATION 

8 

1 

Al.PHA 

0 • 12510-01  +  0,40200-04 

* 

EL*DN 

BETA 

0 » t994D+Ot  +-0,24010-01 

* 

UN 

RMS 

0.34681 17E-01 

8 

2 

ALPHA 

-0,70940-01  +  0,20760-01 

* 

RH 

BETA 

0,13840+01  +-0  >  3036D-01 

* 

RH 

RMS 

0 . 651  5279E-01 

8 

3 

ALPHA 

-0.68330-01  +  0.24950-01 

* 

RH 

+  -0, 

35920+00  *  MP*UA 

BETA 

0. 17370+01  +-0,10080+00 

* 

RH 

+  -0 » ' 

?9<*50-03  *  WNYDN 

RMS 

0.561 1358E-01 

8 

4 

ALPHA 

-0.10180-01  +  0.75190-02 

* 

RH 

BETA 

0,23580+01  +-0 , 1 7770+00 

* 

RH 

RMS 

0.3283286E-01 

8 

5 

ALPHA 

0,12660-02  +  0.78620-03 

* 

RHYER 

BETA 

0.24530+01  +-0,99100-01 

* 

RHYER 

RMS 

0 ,201 1605E-01 

8 

6 

ALPHA 

0.97660-03  +  0,43920-03 

* 

RHYFR 

BETA 

0,22830+01 

RMS 

0.1480237E-01 

8 

7 

ALPHA 

0,25530-02  +-0, 14170-03 

* 

OH 

BETA 

0,22450+01 

RMS 

0 . 21 97240E-01 

8 

8 

ALPHA 

0.60650-01  +-0,92210-07 

t 

RH 

BETA 

0,11070+01 

RMS 

0 « 644531 5E-01 

+-0, 38890+00  *  MPYWA 
+  -0  » 38770-02  *  DNtER 

+  0,55410-01  *  ELtUA 
+  0,35410-01  #  ER 

+  0,14070-01  *  EL*WA 


+  0  •  14430-03  *  DN*F.R 


+  0 1  1 2080-0?  *  El.YDE 


Exhibit  4.3  (Continued) 
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H 

0  H 
N  0 


T 

H 

11 

R 

DF.P 

MAR 

EQUATION 

9 

t 

ALPHA 

-0 . 40240-02 

+  0 . 76890-02 

* 

RH 

BETA 

0.14400+01 

+-0 i 35260-01 

* 

RH4ER 

RMS 

0.6334446E 

-01 

9 

2 

ALPHA 

-0. 4043H+00 

+  0  <  76850-01 

* 

RH 

BETA 

0,  t914D  +  01 

+-0, 11880+00 

* 

RH 

+-0. 19250+00 

* 

ER 

RMS 

0.7961287E 

-01 

9 

3 

ALPHA 

-0.53350+00 

+  0 . 1 0800+00 

* 

RH 

+  -0. 27080+00 

* 

EL*MP 

BETA 

0  >  20180+01 

+-0 . 17350+00 

* 

RH 

+  0  *  t228D-01 

* 

ER 

RMS 

0.5935084E- 

-01 

9 

4 

ALPHA 

-0 .61920-01 

+  0  <  2734D-01 

* 

RH 

+  0.30960-02 

* 

WN*WA 

BETA 

0 . 17290+01 

+-0, 11490+00 

* 

RH 

RMS 

0.4147219E 

-01 

9 

5 

ALPHA 

0,48580-02 

+  0  <  2807D-01 

* 

MP*WN 

+  0.11860-02 

* 

WN*WA 

BETA 

0.19040+01 

+-0,19380-01 

* 

WN*MA 

+-0. 11250-02 

* 

WN*ON 

RMS 

0.2505645E 

-01 

9 

6 

ALPHA 

0.40070-02 

+  0 ♦ 10640-01 

* 

MP*WW 

+  0.94210-03 

* 

WAiON 

BETA 

0 . 1916D+01 

RMS 

0.1763006E 

-01 

9 

7 

ALPHA 

0 . 46790-02 

+  0,35030-01 

* 

MP*WN 

BETA 

0,1 8860+01 

RMS 

0 . 2593861E 

-01 

9 

8 

ALPHA 

0 . 1012D-01 

BETA 

o.?23r.n+oi 

+-0. 13630-01 

* 

HN 

+-0. 10760-02 

* 

RH*WN 

RMS 

0 . 4092572F. 

-01 

.Exhibit  4.3  (Continued) 
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N  0 

T  U  DFP 

H  R  MAR  EQUATION 

10  l  ALPHA  0 i 13290+00 

BETA  0.92590+00 
RMS  0  ,6887224001 

10  2  ALPHA  -0.35690+00  +  0.79740-01  *  RH 

BETA  0.73340+00 

RMS  0.9098264E-01 

10  3  ALPHA  -0.16270+00  +  0.66050-01  *  RH  +-0. 42610+00  *  RH*MP 

+  0,14150+00  *  EL*WA  +  0.33600-02  #  EL*DE 

BETA  0.67450+00  +-0. 2 1780+00  *  FL*WA 

RMS  0  » 4880673E-0 1 

10  4  ALPHA  -0.27830+00  +  0.76420-01  t  RH  +  0,17520-01  *  RH*WA 

+  0,81430-03  *  MN*ER 

BETA  0.23440+01  +-0. 23990+00  *  RH  +-0, 251 70-01  *  RH*WA 

+  -0 ,  t888D-02  *  El.ltWN  +-0, 54590-02  *  ONftDE 

RMS  0 • 3688542E-01 

10  5  ALPHA  0.17930-02  +  0,70480-02  *  RH  +  0.16210-03  *  RH*WN 

+  0,18410-01  *  RH*WA 

BETA  0.15030+01  +-0, 10600-02  *  RH*WH  +-0, 94050-01  *  RH*WA 

+-0,10040-01  *  ON*ER 

RMS  0 . 3120006E-01 

10  6  ALPHA  -0,28950-01  +  0,14500-01  »  RH  +  0,13830-01  *  RH*WA 

BETA  0,19170+01  +-0, 82250-01  *  RH  +-0. 12280+00  *  RH*WA 

+-0, 55620+00  *  MP*WN 

RMS  0.3312006E-01 

10  7  ALPHA  0.58110-02  +  0,11210-01  *  RH 

BETA  -0.18420+00  +  0.248R0+00  *  RH  +  0,41380-01  *  RH*WA 

+  0,11240+00  *  MP*WN  +-0, 28610+01  *  MP*DE 

RMS  0 , 7974254E-0 1 

10  8  ALPHA  0.63710-01 

BETA  0,12300+01 

RMS  0 • 6072881 E-01 


Exhibit  4.3  (Continued) 
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M 

0  H 
N  0 


T 

H 

U 

R 

OF.P 

VAR 

EQUATION 

11 

1 

ALPHA 
+  0. 

-0,80530-01  +  0.25260-01 
70420+00  *  RH*NP 

* 

RH 

+-0 . 6530D-03 

* 

UN 

BETA 

0 . 13400+01  +-0. 21620+02 

* 

MP 

+  -0. 77050-02 

* 

RHT-FL 

RMS 

0.4733886E-01 

11 

2 

ALPHA 

-0.42430+00  +  0 < 7751 0-01 

* 

RH 

+  0.23800+01 

* 

KPtER 

BETA 

0,11200+01  +-0, 28060-02 

* 

EL*WN 

+-0. 22210+00 

* 

MPtDN 

RMS 

0.61  63124E-01 

11 

3 

ALPHA 

0. 16220+00  +  0,61300-01 

* 

El.tWA 

+  0.19040+01 

* 

MPtFR 

BETA 
+  -0. 

0.98150+00  +-0. 46540-03 
30170-02  *  WNtWA 

* 

ELtUN 

+-0. 10560+01 

* 

NPtON 

RMS 

0 » 5082504E-0 1 

11 

4 

ALPHA 
+  0. 

-0.25080+00  +  0,62060-01 
13530+01  *  NP*ER 

t 

RH 

+  0.75980-01 

* 

WA 

BETA 
+  -0, 

0.22130+01  +-0. 18410+00  *  RH  +-0. 16450+00 

82590-03  *  El.tUN  +-0. 68870+00  *  MPHON 

t 

UA 

RMS 

0 , 4053295E-01 

11 

5 

ALPHA 

0.81130-01  +  0,22140-01 

t 

MP*WN 

+  0,25310-0? 

* 

UN*WA 

BETA 
+  -0 , 

0.12140  +  01  +-0,12290-01  *  RH  +-0, 66080-01 

64280-03  *  ELtUN  +-0,56050+00  *  MP*ON 

* 

UA 

RMS 

0 . 3738825E-01 

11 

6 

ALPHA 
+  0, 

0.66700-01  +  0,26340-03 
30590-02  *  WNtUA 

* 

El.tUH 

+  0,17240+00 

* 

MPtDE 

BETA 
+  0. 

0.12060+01  +-0,29500-01 
34350-03  *  WN*ER 

* 

RH*UA 

+-0. 23250-0? 

* 

Ft.tWN 

RMS 

0 . 3872498E-01 

11 

7 

ALPHA 
+  0, 

-0,24530+00  +  0,52490-01 
15220-02  *  WN*UA 

* 

RH 

+  0 , 97990+00 

* 

MPtER 

BETA 

0.12030+01  +-0. 41470-01 

* 

FL.tER 

RMS 

0.5391571E-01 

11 

8 

ALPHA 

0.10580+00 

BETA 

0.11360+01  +-0. 15680-02 

t 

ELtUN 

RMS  0.5641866E-01 


Exhibit  4.3  (Continued) 
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M 

0  H 

N  0 

T  II  PEP 

H  R  OAR  EQUATION 

12  1  ALPHA  0 » 1087D+00  +  0.41240-01  *  El.  +-0,19230+00  t  RH*HP 

BETA  0  *  11250+01  +-0.A103D-01  *  ELCER  +  0.153tD+01  *  MP*DE 

RMS  0.58751 15E-01 

12  2  ALPHA  0.1 5180+00  +  0,49770-03  *  EL.CWN  +-0 , 244811-02  »  ELCPE 

BETA  0 . 10A30+01  +-0.1822)1-01  *  ONCER 

RMS  0 » A859B81 E-01 

12  3  ALPHA  0 . 1 3830+00  +  0.24A3D-03  *  ELCWN  +  0,22530-03  *  UNION 

BETA  0 . 11310+01  +-0.34A2D-01  «  DNCER 

RMS  0.4941 199F-01 

12  4  ALPHA  0.1A25D+00  +  0.25A4D-03  *  UNCDN 

BETA  O.JUOn  +  Ol  +-0 . 30830-01  *  ONCER 

RMS  0  >  48393A8E-01 

12  5  ALPHA  0 . 1 1 19P+00  +  0.14340-02  *  WNCWA  +  0,21930-03  C  UNCPN 

BETA  0.11810+01  +-0. 10490-02  *  UNCDN 

RMS  0 . 4A28299E-01 

12  A  ALPHA  0,10720  +  00  +  0.1387D-03  *  EI.CWN  +  0.28A30+00  C  ELCUA 

+  0,19840-03  *  UNCDN 

BETA  0,12030  +  01  +-0.8A17D+00  *  ELCWA  +  0,30150-01  *  ONCER 

RMS  0.42AA235E-01 

12  7  ALPHA  0.10830  +  00  +  0.47570-03  *  RHCON  +  0.31310-03  *  El.CUN 

BETA  0,12090+01  +-0.29AAD-01  *  ONCER 

RMS  0.5020483E-01 

12  8  ALPHA  0,10970+00  +  0.1A420-01  C  EI.CER 

BETA  0.11700+01  +-0, 79120-03  f  UNCDN 

RMS  0 i 7029019E-01 
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Exhibit  4.4 


Exhibit  4.5 


Exhibits  4.4  and  4.5  indicate  the  month  and  hour  period  for  which  the 
single  factor  or  interaction  terms  are  used  in  the  prediction  equations  for 
a  and  0  respectively.  In  Exhibits  4.4  and  4.5  the  variables  are  referred  to 
by  index  numbers  1  to  36  for  economies  of  space.  Exhibit  4.6  gives  the  coded 
terms  corresponding  to  these  index  numbers. 
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Exhibit  4.6 

TERMS  CORRESPONDING  TO  INDEX  VARIABLES  IN  EXHIBITS  4.4  AND  4.5 
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Exhibit  4.7 

FREQUENCY  OF  APPEARANCE  OF  VARIABLES  IN  EXPRESSIONS  FOR 
ALPHA  AND  BETA  ("VARIABLES  MODEL") 


Exhibit  4.7  gives  the  frequency  of  appearance  of  variables  in 
expressions  for  a  and  0.  It  is  worth  remarking  that  relative  humidity  appears 
more  often  in  expressions  for  predicting  a  and  6  than  any  other  variable  used, 
and  as  a  "main  effect"  term  it  appears  in  29  of  the  96  expressions  for  a  and 
16  of  the  96  expressions  for  0.  Elevation  and  wind  also  appear  in  many  of  the 
equations,  but  almost  always  as  "interaction”  terms.  This  suggests  that 
humidity  by  itself  is  useful  in  obtaining  expressions  for  a  and  0,  but 
elevation  and  wind  are  useful  only  in  combination  with  other  variables. 

Interaction  terms  which  appear  frequently  in  expressions  for  a  are 
EL  x  VN  and  EL  x  WA,  each  of  which  occur  in  15  of  the  expressions.  Those 
which  appear  frequently  in  expressions  for  0  are  EL  x  WN,  WN  x  DN,  EL  x  WA, 

RH  x  WA,  RH  x  WN,  WN  x  WA,  and  DN  x  ER  which  occur  41,  36,  22,  20,  15,  15  and 
15  times  respectively.  The  interaction  terms  WN  x  DN  and  WN  x  EL  occur  most 
frequently  in  the  expressions  for  8  for  the  Winter  and  Spring  months. 

There  are  situations  where  no  knowledge  of  observable  climatological, 
geographical  or  other  variables  is  known  (or  we  wish  to  ignore  it  in  the 
interest  of  simplicity).  For  these  cases,  we  have  developed  constants  models 
for  the  entire  region  of  West  Germany  using  the  Weibull  distribution,  which 
have  fixed  values  for  a  and  8  for  each  of  the  96  combinations  of  month  and 
three-hour  period.  Exhibit  4.8  gives  values  for  the  constants  and  also  RMS 
values  over  all  stations  and  distances.  These  values  were  obtained  by  the 
same  method  as  for  the  previous  models. 
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Exhibit  4.8 

VALUES  OF  o,  6  AND  RMS  FOR  "CONSTANT  MODEL" 
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The  expression  minimized  was  (3.1)  of  section  3.0,  with  expressions  (3.3) 
and  (3.4)  containing  only  the  constant  terms  Yq  and  Sq  respectively.  In 
other  words,  the  best  (a, 6)  combination  for  the  entire  region  is  found  by 
minimizing 

30  14  r  -i  2 

S  E  E  (x  )  -  F(x  ;o,B)  (4.1) 

j-1  i=l  L  2  1  J 

with  respect  to  a  and  8.  (x^)  is  the  empirical  probability  of  less  than  x^ 

miles  of  visibility  at  station  j.  F(x^;o,8)  is  the  Weibull  model  for 
visibility.  Note  that  in  this  case  we  consider  (o,8)  as  constants  for  the 
entire  region  of  West  Germany.  That  is,  we  obtain  a  single  regional 
probability  model  for  a  given  month  and  hour  period. 

For  example,  for  February  for  the  third  hour  period  we  have  the 
regional  model  given  by 

F(x)  *  1  -  e  (4.2) 

This  regional  model  can  be  used  for  the  entire  West  German  region  for 
February  third  hour  period  and  yields  an  overall  RMS  of  .078. 


5.0  ACCURACY  OF  THE  RESULTS 

Exhibit  5.1  gives  a  comparison  of  the  accuracies  of  the  two  models. 
The  values  given  are  the  RMS  values  of  the  differences  between  the  model 
predicted  probabilities  and  the  observed  probabilities  (RUSSWO  values)  over 
all  distances  and  over  all  of  the  thirty  stations. 
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Exhibit  5.1 

RMS  VALUES  FOR  THE  "VARIABLES  MODEL"  AND  THE  "CONSTANTS  MODEL" 
(The  bottom  figure  for  each  pair  is  for  the  Constants  Model) 
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It  is  worth  noting  that  the  RMS  values  are  lowest  for  the  month  of  May 
and  for  hour  period  6  (1500-1700  hours) .  The  RMS  values  are  highest  for 
January  and  for  hour  period  2  (0300-0500  hours) .  This  may  be  because 
visibilities  are  worst  for  these  times,  or  for  the  case  of  hour  period  2, 
that  the  observations  are  less  accurate. 

There  are  some  apparent  anomolies  in  Exhibit  5.1.  For  March,  and  the 
first  three  hour  period,  the  RMS  value  for  the  Variables  Model  is  larger 
than  the  RMS  value  for  the  Constants  Model.  The  explanation  is  that  in  the 
development  of  the  Variables  Model,  for  a  specific  hour  period  some  of  the 
stations  had  to  be  omitted  because  values  of  a  climatological  variable  were 
not  available  for  that  station.  For  the  period  0000-0200  hours  for  March, 
the  stepwise  procedure  called  for  mean  monthly  windspeed  and  mean  monthly 
precipitation.  For  two  of  the  21  stations  with  RUSSWO  data,  the  latter  was 
not  available  and  for  another  the  former  could  not  be  obtained.  These 
stations  were  ones  for  which  the  Constants  Model  gave  better  than  average 
fits. 

Exhibit  5.2  compares  the  RMS  for  the  30  German  stations  for  January 
1000  hours  for  the  Constants  Model,  the  Variables  Model  and  the  individual 
station  models.  The  overall  RMS  values  for  the  three  models  are  .099, 

.074,  and  .031.  Again  RMS  values  are  not  given  for  the  variables  model  for 
a  number  of  stations  because  values  for  the  variables  called  for  by  the 
stepwise  regression  model  were  not  available.  Availability  of  these  values 
would  be  expected  to  improve  the  RMS  value  of  .074.  It  may  be  noted  that  the 
RUSSWO  periods  of  observation  for  Baumholder  and  Siegenberg  are  6  years  and 
3  years  respectively,  which  may  account,  at  least  partially  for  their  large 


RMS  values. 


29 


RMS  VALUES 


CONSTANTS 

VARIABLES 

INDIVIDUAL 

STATION 

MODEL 

MODEL 

STATION  MODELS 

Hahn 

.099 

.060 

.013 

Bitburg 

.052 

.046 

.014 

Ramstein 

.107 

.090 

.012 

Spangdahlem 

.061 

.045 

.011 

Tempelhof 

.078 

.078 

.030 

Ansbach 

.028 

.026 

.022 

Fulda 

.026 

.021 

.020 

Erding 

.082 

.076 

.009 

Feucht 

.041 

.047 

.021 

Baumholder 

.209 

.215 

.051 

Bad  Kreuznach 

.023 

— 

.012 

Bad  Tolz 

.171 

.132 

.032 

Zweibrucken 

.040 

.034 

.017 

Wiesbaden 

.022 

.046 

.010 

Finthen 

.051 

.036 

.011 

Furth 

.133 

— 

.048 

Hanau 

.035 

— 

.025 

Gablingen 

.053 

.062 

.033 

Giebelstadt 

.066 

.051 

.014 

Grafenwohr 

.077 

.085 

.022 

Heidelberg 

.045 

.020 

.012 

Illeshein 

.086 

— 

.016 

Kit z ingen 

.060 

.032 

.009 

Nurnberg 

.071 

— 

.020 

Coleman 

.069 

— 

.016 

Wertheim 

.121 

.046 

.016 

Schwaebisch  Hall 

.082 

.067 

.033 

Sembach 

.026 

.025 

.013 

Siegenberg 

.296 

— 

.116 

Echterdingen 

.072 

.059 

.021 

Overall 

.099 

.074 

.031 

Exhibit  5.2 


RMS  OF  FITS  FOR  WEST  GERMAN  STATIONS 
January  1000  Hours 


It  is  worth  noting  that  the  RMS  values  are  differences  between  the 
model  estimates  of  the  probabilities  of  visibility  and  those  of  the  observations 
from  the  RUSSWO  data.  The  latter  are  based  on  the  observers  subjective 
classification  which  contribute  to  an  inherent  error  which  in  turn  contributes 
to  the  increased  size  of  all  RMS's.  It  is  indeed  quite  possible  that 
estimates  from  some  of  the  models  (for  example  those  of  the  individual 
stations)  might  be  an  improvement  over  the  corresponding  RUSSWO  estimates. 


6.0  USING  THE  VISIBILITY  MODELS 


Example  1  It  is  desired  to  obtain  the  probability  of  visibility  less 
than  A  miles  at  a  location  where  there  are  no  historical  records  for  visibility. 
The  probability  is  needed  for  November  at  approximately  0700  hours.  We  are 
willing  to  estimate  the  following  for  the  location  (it  is  not  near  a  major 
body  of  water) . 

November  mean  precipitation  in  inches  1.5 

Elevation  in  feet  1000 

Mean  relative  humidity  (per  cent) 

for  0700  hours  in  November  80 

North-South  elevation  differential  (feet)  -500 

East-West  elevation  differential  (feet)  100 

Mean  November  windspeed  (knots)  A 

Relative  elevation  (elevation/average  elevation)  1.1 

From  Exhibit  A. 3,  we  estimate  o  and  g  for  the  problem  using 
a  -  .1622  +  .0613  EL*WA  +  1.90A  MP*ER 


3  -  .9815  -  .OOOA65A  EL*WN  -  1.056  MP*DN  -  .003017  WN*WA 
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U^irig  Exhibit  4.2, 

EL  -  1000 3/109 
WA 

MP  -  1.53/1000 
ER  =  l.l2 
WN  -  43/10 
DN  =  (-500)2/l02 


=  1 
=  0 

=  .003375 
=  1.21 
=  6.4 
=  2.5 


Calculating,  we  have  a  =  .1700 

8  =  .9696 


Using  the  Weibull  model. 


Prob  f  Visibility  <  4 


]  =  !  -  e"1700(4> 


.9696 


.48 


The  probability  of  visibility  less  than  4  miles  is  estimated  as  .48. 


Example  2  We  wish  to  find  the  probability  of  visibility  less  than 
1  mile  in  March  at  1000  hours  for  a  specified  location.  There  are  no 
historical  records  for  visibility.  We  choose  not  to  use  estimates  for 
geographic  or  climatic  variables.  Using  Exhibit  4.8  we  find  that  for 
March  and  the  fourth  3-hour  period,  the  values  for  a  and  8  can  be  taken  as 
.078  and  1.32  respectively.  The  probability  of  visibility  less  than  or 
equal  to  x  miles  is  given  as 


1 


-,078x 


1.32 


e 


The  probability  that  visibility  is  less  than  1  mile  is  estimated  to  be  .075. 
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7.0  ADDITIONAL  COMMENTS 

In  the  present  report  the  RMS  values  (estimates  of  error  for  calculated 
probabilities)  are  internal  estimates.  A  better  estimate  of  error  could  be 
obtained  using  data  from  a  set  of  stations  not  used  in  determining  the 
prediction  equations.  Future  work  includes  use  of  a  new  set  of  data  from 
several  stations  in  West  Germany  to  evaluate  the  capability  of  the  model  for 
these  stations. 

In  addition,  we  plan  to  utilize  sample  re-use  techniques  such  as  those 
described  by  Geisser  (1975). 

It  will  be  necessary  to  develop  models  using  data  from  other  regions.  It 
is  possible  that  a  small  set  of  models  might  be  sufficient  for  a  major  portion 
of  the  globe.  The  present  work  should  be  extended  to  other  weather  elements 
such  as  ceiling,  windspeed  and  precipitation. 


d 
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APPENDIX 


STATIONS  USED  IN  THE  STUDY 


NO. 

STATION 

WBAN 

LAT. 

LONG. 

ELEV.  (FEET) 

1 

Hahn  AB 

34055 

49.95 

7:27 

1650 

2 

Bltburg  AB 

34049 

49.95 

6.57 

1228 

3 

Ramstein  AB 

34050 

49.43 

7.58 

780 

4 

Spangdahlem  AB 

34054 

49.97 

6.70 

1196 

5 

Tempelhof  APRT 

35104 

52.47 

13.40 

164 

6 

Ansbach  AAF 

34172 

49.32 

10.63 

1542 

7 

Fulda  AAF 

35053 

50.53 

9.63 

1010 

8 

Erdlng  AS 

34168 

48.32 

11.93 

1522 

9 

Feucht  AAF 

34198 

49.38 

11.18 

1265 

10 

Baumholder  AAF 

34077 

49.65 

7.30 

1408 

11 

Bad  Kreuznach  AAF 

34070 

49.87 

7.88 

355 

12 

Bad  Tolz  AAF 

34197 

47.77 

11.60 

2360 

13 

Zweibrucken  AB 

34058 

49.22 

7.40 

1132 

14 

Wiesbaden  AB 

35010 

50.05 

8.33 

470 

15 

Finthen  AAF 

34075 

49.97 

8.15 

769 

16 

Furth  AAF 

34176 

49.50 

10.93 

1000 

17 

Hanau  AAF 

35009 

50.17 

8.95 

377 

18 

Gabl ingen  AAF 

34196 

48.45 

10.87 

1530 

19 

Giebelstadt  AUX  AF 

34036 

49.67 

9.88 

985 

20 

Grafenwohr  AAF 

34189 

49.70 

11.95 

1370 

21 

Heidelberg  AAF 

34046 

49.40 

8.65 

369 

22 

Illeshelm  AAF 

34190 

49.47 

10.38 

1060 

23 

Kitzinger  AAF 

34191 

49.75 

10.20 

699 

24 

Nurnberg 

34177 

49.50 

11.08 

1053 

25 

Coleman  AAF 

34068 

49.57 

8.47 

334 

26 

Wertheim  AAF 

34076 

49.77 

9.48 

1120 

27 

Schwaebisch  Hall  AAF 

34074 

49.17 

9.78 

1303 

28 

Sembach  AB 

34056 

49.52 

7.87 

1052 

29 

Siegenberg  Gunnery  Range 

34199 

48.75 

11.80 

1325 

30 

Echterdingen  ARPT 

34041 

48.68 

9.22 

1306 

